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Abstract 



We investigate the kinetics of many-species systems with aggregation of 
similar species clusters and annihilation of opposite species clusters. We find 
that the interplay between aggregation and annihilation leads to rich kinetic 
behaviors and unusual conservation laws. On the mean-field level, an exact 
solution for the cluster-mass distribution is obtained. Asymptotically, this so- 
lution exhibits a novel scaling form if the initial species densities are the same 
while in the general case of unequal densities the process approaches single 
species aggregation. The theoretical predictions are compared with numerical 
simulations in ID, 2D, and 3D. Nontrivial growth exponents characterize the 
mass distribution in one dimension. 
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Irreversible aggregation and annihilation processes occur in many natural phenomena 
]1].|2|]. The kinetics of each of these processes is well understood both on the mean- field 
level which provides an accurate description when fluctuations in reactant densities can be 
ignored and in the fluctuation-dominated regime which occurs in low-dimensional 

systems, see ||, and references therein. It was found in particular that the cluster- mass 
distribution in aggregating systems typically approaches a universal scaling form [§,§J . 

In this work we investigate the competition between aggregation and annihilation, which 
gives rise to a surprisingly rich kinetic behavior. Our model is the simplest implementation 
of the two underlying processes, aggregation and annihilation: Similar species clusters com- 
bine to form larger clusters while dissimilar clusters combine to form an inert substitute. It 
is worth noting that when the number of species n diverges, n — > oo, our model is equivalent 
to single-species annihilation, while in the other extreme, n = 1, single-species aggrega- 
tion is recovered. Hence, the process is well suited for investigating the interplay between 
aggregation and annihilation. We find that the exponents characterizing the kinetics are 
nonuniversal in that they may depend on the reaction rates as well as the initial conditions. 
Still, there is a basic scaling form which describes the process in the long-time limit. 

We shall study our model both in the mean-field limit and in the diffusion-controlled 
limit for d = 1, 2 and 3. The mean-field approach to the binary reaction process assumes 
that the reaction proceeds with a rate proportional to the product of the reactants densi- 
ties. Thus the mean-field approximation neglects spatial correlations and therefore typically 
holds in dimensions larger than some critical dimension d c . For n-species pure annihilation 
processes, it was suggested || that d c = A{n — l)/{2n — 3), the result which turns into rigor- 
ously established values of the critical dimension for two- and single-species annihilation, 
d c — 4 and 2, respectively (the latter case corresponds to n — oo). For pure aggregation 
processes, the critical dimension depends on the details of the reaction events and on the 
relation between the diffusion coefficient and the mass of the cluster ||. Numerically, we 
investigate the particle coalescence model (PCM) in which clusters occupy single lattice sites 
and perform nearest neighbor hoping with the diffusion coefficient independent of the mass. 
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For the PCM, it is known that d c = 2 ||. Since the aggregation-annihilation model inter- 
polates between single-species aggregation and single-species annihilation processes, both of 
which belong to the same universality class, it is natural to expect that the critical dimen- 
sion remains the same, d c = 2. We confirm the mean-field predictions above this critical 
dimension numerically. In one dimension, spatial correlations are relevant asymptotically. 
A simple heuristic argument provides the exact asymptotic behavior of the concentration 
and a good approximation for the growth of the typical mass. 

Following the above discussion, similar species aggregation is described by the binary 
reaction scheme 

Ai + Aj — > A,i + j, Bi + Bj — > Bi + j, (1) 

where Ai denotes a cluster consisting of i monomers of species A, and similarly for B^. 
The two-species case contains the generic many-species behavior. Therefore, we shall focus 
on the two-species situation and cite the many-species results if appropriate. Annihilation 
between dissimilar clusters reads 

Ai + Bj -> inert. (2) 

Thus, we assume that dissimilar clusters completely annihilate independent of their masses. 
In some situations it may be more reasonable to assume a partial annihilation, where the 



monomer difference number i — j is conserved [jrU-12 . 

Denote by ak{t) and bk(t) the densities of A- and -B-clusters, consisting of k monomers. 
Then the mean-field rate equations for the two-species aggregation-annihilation process read 

«fc = H a i a j ~ 2a k(a + b), b\= bibj -2b k (a + b). (3) 

i+j=k i+j=k 

Here a(t) and b(t) denote the total densities of A- and -B-clusters, a = J2k>i a k an d b = 
Sfc>i bk, and the overdot denotes the time derivative. In writing Eq. (3) we have assumed 
that the aggregation rate equals the annihilation rate. We set this constant rate to unity 
without loss of generality. 
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Summing up Eq. (3) we obtain 

a = -a 2 - 2ab, b = -b 2 - 2ab. (4) 

For the symmetric initial conditions, a(0) = 6(0), we get a = —3a 2 which is immediately 
solved to find 

For the asymmetric initial conditions, a(0) > 6(0), it is helpful to rewrite Eq. (4) in terms 
of u = (a + b)/(a — b) and v = a — b. This gives 

u = — v (u 2 — 1), v — —uv 2 . (6) 
2 

Expressing v as a function of u yields dv/du = 2uv/ {u 2 — 1), and as a result 

1 >' 4a(0)6(0) ' ( ' 

Inserting Eq. (7) into Eq. (6), we obtain a closed equation for u = u(t) with the solution 



2u 2u(0) 

+ ln 



u 2 - 1 u 2 (0) - 1 



u - 1 u(0) + 1 
u+ 1 m(0) - 1 



2ift. (8) 



The densities a(t) and 6(t) of the total number of A- and _B-clusters, a = v{u + l)/2 and 
6 = f (u— l)/2, cannot be expressed as explicit functions oft. However, asymptotically they 
exhibit simple power-law behaviors 

a~r\ 6~a(0)6(0)[a(0)-6(0)]~ 3 r 2 for t -> oo. (9) 

This indicates that the majority and the minority species evolve very differently. Similarly, 
in the many-species case the total density of the majority species decays as t^ 1 while the 
minorities densities decay as t~ 2 . 

We turn now to the determination of the cluster densities. Introducing the generating 
functions, 

oo oo 

A(z,t) = J2z'a J (t), B(z,t) = J2^b J (t), (10) 



transforms the governing equations into 



A = A 2 -2A(a + b), B = B 2 -2B(a + b). 



(11) 



These Bernoulli equations are straightforwardly solved to get 



A(z,t) = 



A (z)E(t) 



l-A (z)f*dt>E(t'y B{Z,t) l-B (z)J*dt>E(t>Y 

with the shorthand notations A (z) = A(z,t = 0), B (z) = B(z, t = 0), and E(t) = 
exp (-2 f* dt'[a(t') + b(t')}). 

Eq. (12) represents the general solution for arbitrary initial conditions. Consider now 
the simplest but important case of monodisperse, generally asymmetric, initial conditions 



B (z)E(t) 



(12) 



afc(o) = 4i, h(0) = \S kl . 



(13) 



These initial conditions imply A (z) = z, B (z) = \z. Expansion of the resulting generating 
functions A(z,t) = zE{t)[l- z ^dt' E{t')}- 1 and B(z,t) = XzE(t)[l- zX f*dt' 'E{t')]~ l yields 

k-i 



a k (t) = E{t) f dt'Eit') 



k 1 , b k (t) = X k E(t) 



( dt'Eit') 
Jo 



(14) 



In the symmetric case A = 1 we easily compute E(t) — (1 + 3t) 4//3 and as a result 



ik-i 



a k (t) = b k (t) = (1 + 3t)- 4 / 3 [l - (1 + Sty 1 / 3 



Asymptotically, the cluster-mass distribution approaches the scaling form 



(15) 



a k 



(t) ~r 4 / 3 exp(-:r), x ~ kt~ l l\ 



(16) 



The total mass densities, m a (t) = J2k>i ka k (t) and rribit) = J2k>i kb k (t), decrease with time, 
= m b ~ t~ 2 / 3 . Note also that the quantity / = m a a~ 2 / 3 is conserved by the dynamics of 
the aggregation-annihilation model subject to arbitrary symmetric initial conditions. The 
conservation of I is verified by a direct computation which makes use of the evolution 
equations, d = —3a 2 and m a = —2am a . For pure aggregation processes, the total mass is 
conserved. Thus the quantity / plays the role of a "mass" in the aggregation-annihilation 



model. Unfortunately, we could not understand the physical meaning of the conservation 
law I = const. 

In the asymmetric case we shall set A < 1 thus forcing 5-species to be a minority. Making 
use of Eq. (3) one can express E(t) = exp (-2J*dt'[a(t') + 6(f)]) as an explicit function of 
a and b, E = (a(0) _1 — 6(0) _1 )/(a _1 — b^ 1 ). In the long-time limit we use the asymptotic 
values, a ~ t~ l and b ~ A(l - A)~ 3 t -2 , to compute E{t) ~ A _1 (l - \)b ~ (1 - A)~ 2 r 2 . 
Therefore, 

a fc (t)~(l-A)-V 2 exp(-y), b k (t) = \ k a k (t), y= (17) 

Thus, for arbitrary initial conditions the majority species cluster-mass distribution can be 
written in the scaling form 

a k (t)~r w $[k/S(t)}, S(t)~t z , (18) 

where S(t) is the characteristic mass. The scaling function is exponential, $(a;) = exp(— x), 
both for symmetric and asymmetric initial conditions. In the symmetric case, the governing 
exponents are w = 4/3 and z — 1/3. The asymmetric case is equivalent to single species 
aggregation, and the exponents are w = 2 and z — 1. The minority, on the other hand, 
does not scale according to the usual definition although it can be expressed in the modified 
scaling form, b k (t) ~ A k t~ w &[k/ S(t)}. The modified scaling form also indicates that two 
different mass scales are associated with the minority species. A growing scales S(t) ~ t, 
which is forced by the majority species, and a time independent scale S\ — 1/(1 — A) which 
dominates in the long-time limit. The latter scale diverges, S\ = 1/(1 — A), in the limit 
A — > 1. It is also instructive to compute the total mass densities. By summing Eq. (17) we 
find that as t — > oo, 

m a (t) -> (1 - A) 2 , m b (t) -> A(l - XyH' 2 . (19) 

Thus the final mass difference Am^ = m a (oo) — mfe(oo) = m a (oo) may be expressed through 
the initial mass difference Am = m a (0) — m&(0) = 1 — A via a surprisingly simple relation, 
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AttIoo = (Am ) 2 . In comparision, for the pure annihilation process, A4 + Bj — > 0, the 
final mass density is significantly larger, Am^ = Am . More generally, the concentration 
difference, c_ = a — b, is a conserved variable in the annihilation process, as can be seen from 
the evolution equations a = b = —lab. In the aggregation-annihilation process a related 
"hidden" conservation law exists. From the rate equations, one can verify that the quantity 
c_ = m a mi>C-/ab is conserved. This unusual conservation law is trivially satisfied in the 
case of equal initial densities, where the aforementioned aggregation-like conservation law 
m a a~ 2 / 3 = const holds. We conclude that the aggregation-annihilation process exhibits two 
nontrivial conservation laws which are generalization to the usual conservation laws that 
underly aggregation and annihilation separately. 

We turn now to the general n-species model. For asymmetric initial conditions the 
majority species scales as in single-species aggregation, while the minorities scale only in 
a modified sense. For the symmetric initial conditions, the cluster-mass distribution ap- 
proaches the scaling form of Eq. (18) with 

2n , 1 

W n = o 7> aIld Z n = o r ( 20 ) 

2n — 1 2n — 1 

Thus the exponents depend on the number of species n. When n — 1, single-species aggre- 
gations is recovered w\ — 2 and z\ — 1, while the limit n — > 00, corresponds to single-species 
annihilation, w — 1 and z — 0. Denoting by a 1 (m % ) the concentration (mass) of the i th 
species, then the quantity, M{a % — a^)/a l a^ is conserved for every i 7^ j independent of the 
initial conditions. The factor M is given by M n ~ l = Il™ =1 a\ Only n — 1 of these conserva- 
tion laws are independent, similar to the n — 1 independent conservation laws that unedrly 
the pure annihilation process, a % — a- 7 = const, for i ^ j. For the case of symmetric initial 
conditions, the aggregation induced conservation laws read rn l {a l )~^ 2n ~ 2 ^^ 2n ~ 1 ^ = const for 
i — 1, . . . , n. 

Generally, we expect that the exponents characterizing the scaling behavior are universal, 
i. e., they depend only on important aspects of the kinetics. The fact that the exponents do 
depend on the number of species does not contradict universality. However, the exponents in 



the aggregation-annihilation system may depend also on the reaction rates. To demonstrate 
that we consider the general case where the aggregation rate and the annihilation rate are 
different. We set the aggregation rate to unity as previously and denote the annihilation 
rate by J. In the n-species case with symmetric monodisperse initial conditions we have 

a k = a^j - 2a k [l + (n - l)J]a, a k (t = 0) = 8 kl . (21) 

i+j=k 

By employing the above technique we solve Eq. (21) and find 



k-l 



a k (t) = (l + ut)- 1 - 1 ^[l-{l + ut)- 1 ^\ , (22) 

where the notation v = 1 + 2(n — 1) J has been used. The cluster mass distribution obeys 
the scaling form of Eq. (18) with the exponents 

l + (n-l)J , 1 

w n = 2 ; r — and z n = ; — . (23) 

l + 2(n-l)J l + 2(n-l)J V ; 

Interestingly, both exponents dependent continuously on the (relative) magnitude of the 
annihilation rate, J. This nonuniversality contrasts the bulk of previously investigated 
aggregation models ||] (see, however, ]T3|Jn|). Varying the rate ratio represents another 
way of interpolating between aggregation and annihilation. Indeed, independent of n, when 
J = 0, single-species aggregation is recovered, while the limit J —>■ oo corresponds to single- 
species annihilation. 

We now consider aggregation-annihilation in the diffusion-controlled limit in low dimen- 
sions. We shall study the particle coalescence model (PCM) in which clusters occupy single 
lattice sites, hop to nearest neighbor cites with a rate independent of the mass, and aggre- 
gate or annihilate instantaneously whenever they meet. The PCM allows one to focus on the 
kinetic aspects of the process and is well suited for numerical implementation. Additionally, 
in the mean-field approximation the PCM is governed by the same Eq. (5) which has been 
examined previously. 

Let us ignore the mass and the identity of the clusters and denote an arbitrary cluster by 
C . Then the reduced reaction process can be described by the reaction scheme C + C — » C 
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and C + C — > for aggregation and annihilation, respectively. For both processes, the total 
density c = Y^k>\ c k is inversely proportional to N(t), the number of distinct sites visited by 
a single random walk in d dimensions 0. N(t) is well known in probability theory [|14]] and 
thus the density c, c ~ iV _1 , behaves as 

f r 1 / 2 d = 1; 
c(t) ~ | r 1 log(t) d = 2; (24) 

U" 1 d>2. 

In other words, the time variable t is replaced with a modified time variable N(t). We further 
assume that the rate equation theory describes the process in low dimensions, with the time 
variable N(t). Although this approach is a heuristic one, it provides a good approximation 
for the subcritical behavior. For the symmetric initial conditions, the characteristic mass 
for the case n = 2 is given by 

r *V6 d = 1; 

S{t)~\ (t/log(t)) 1/3 d = 2; (25) 
1 1 1 / 3 d > 2. 

Using the asymptotic forms of the concentration and the typical mass, a scaling form similar 
to the one of Eq. (18) can be written. The same analysis can be repeated for the general 
multi-species case, and we quote the scaling exponents in one dimension only, 

w n = , and z n = — -. (26) 

2n-l' 2(2n-l) V ; 

In the extreme cases of n = 1 and n — > oo, these exponents agree with the solutions to 

single-species aggregation and single-species annihilation, respectively Hl~5| . In the case of 



asymmetric initial conditions, the majority species concentration decays according to single- 
species aggregation, or equivalently, according to Eq. (23). The minority species, on the 
other hand, decays much faster. Using the above heuristic argument, we find that in the 
long-time limit, the minority-species concentration decays as t~ l in ID, and as t -2 log 2 (t) in 
2D. 

It is interesting to compare the above theoretical predictions with numerical simulations. 
We have performed simulations for d = 1,2, and 3. The numerical implementation of the 



process is simple. Initially all L d sites on the cubic lattice are occupied with monomers. 
An elemental simulation step consists of picking a cluster at random and moving it to a 
randomly chosen neighboring site. If the site is occupied, an aggregation or an annihilation 
event takes place, depending on the identity of the two clusters. Time is updated by the 
inverse of the total number of particles in the system after each step. The linear dimension 
of the lattice used in the simulation was L = 3 x 10 7 , 4 x 10 3 , 2 x 10 2 in ID, 2D, and 3D 
respectively. The corresponding number of particles was thus 30 x 10 6 , 16 x 10 6 , and 8 x 10 6 
for ID, 2D, and 3D. 

We discuss first results for the symmetric initial conditions involving two species. To 
verify that the process belongs to the PCM universality class we measured the density of 
particles as a function of time. Indeed, Eq. (24) appears to hold asymptotically (see Figure 
1). We also measured the average cluster size (k(t)) = J2k c k/ J2 c k- We expect that the 
average cluster size is proportional to the characteristic size S(t), asymptotically. For d = 3, 
the mean-field prediction, {k(t}) ~ t 1 ^ is verified (see Figure 2). Note that the simulation 
results are reliable up to time oc L 2 due to finite size effects. At the critical dimension, d c = 2, 
the average mass grows slightly slower, consistent with the logarithmic correction of Eq. (25). 
In ID, the growth exponent z as determined by a least square fit is z = 0.190 ± 0.002. 
However, the corresponding value from the approximate theory is lower z = 1/6. The 
simulation was carried over a relatively large temporal range, suggesting that the growth 
exponent is indeed different than suggested by the heuristic theory. A similar trend is 
observed when more than two species are involved. We find that the numerically determined 
exponents z 3 = 0.12 and z^ = 0.09 are slightly higher than their theoretical counterparts 
zs = 1/10 and z A = 1/14. We have made some additional consistency checks. For example, 
we verified that the cluster mass distribution follows the scaling form of Eq. (18). Above 
the critical dimension the scaling function is indeed a simple exponential, in agreement with 
the rate equation predictions. 

In summary, we introduced an annihilation-aggregation model and presented the solution 
to the governing rate equations. The mass distribution follows a general scaling form with 
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nontrivial growth exponents, and the process is characterized by unusual conservation laws. 
The kinetic behavior depends both on the number of species and on the specific rates of 
the annihilation and the aggregation processes. Numerical simulations confirm the rate 
equation predictions above the critical dimension. Below the critical dimension, a heuristic 
argument provides a good estimate for the mass distribution. Nontrivial exponents describe 
the process in one-dimension and it would be interesting to study this system rigorously. 
Furthermore, the underlying conservation laws in low dimensions are especially intriguing. 
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Figure Captions 

Fig. 1 The total density versus time in ID, 2D and 3D. Shown are the simulation data 
(bullets) and the theoretical prediction of Eq. (24) (solid line). 

Fig. 2 The average cluster size (k(t)) versus t in ID, 2D, and 3D. Shown are simulation 
data (bullets) and lines of slope 0.19 (solid line) and 1/3 (dashed line) for reference. 
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